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We present data from Fermilab experiment E781 (SELEX) on the hadroproduction asymmetry 
for A c compared to A+ as a function of if and p\ distributions for A+ . These data were measured 
in the same apparatus using incident 7r~, E~ beams at 600 GeV/c and proton beam at 540 GeV/c. 
The asymmetry is studied as a function of xp. In the forward hemisphere with xf > 0.2 both baryon 
beams exhibit very strong preference for producing charm baryons rather than charm antibaryons, 
while the pion beam asymmetry is much smaller. In this energy regime the results show that beam 
fragments play a major role in the kinematics of A c formation, as suggested by the leading quark 
picture. 
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I. INTRODUCTION 

The conventional leading order (LO) and next to lead- 
ing order (NLO) charm-anticharm production diagrams 
at the quark level show little or no asymmetry in the xf 
or p t behaviour of the quark and antiquark, as viewed in 
the laboratory frame. Previous studies, primarily using 
pion beams at a variety of energies up to 500 GeV/c, 
have shown that in many cases charm mesons that share 
a quark in common with the beam hadron, so-called 
leading charm, are produced more copiously in the for- 
ward hemisphere than those states which do not have a 
shared quark, i.e., non-leading charm M. Recently, this 
has been extended to charm baryon production by a £~ 
beam, where asymmetry for A+ compared to A c as a 
function of xf was reported to be large There are 
limited-statistics proton data available for central pro- 
duction, where only a 90% confidence level lower limit is 
reported 

Because the quark-level processes are thought to be 
symmetric, the asymmetries are viewed as features of the 
hadronization process Q or as a manifestation of an in- 
trinsic charm content of the beam hadron J5| . In this ex- 
periment the A+ shares a quark in common with all three 
beam hadrons, while the antibaryon has an antiquark in 
common with the 7r~ beam only. We shall compare the 
hadroproduction characteristics for these charm baryon 
and antibaryon states to illuminate the role of the beam 
fragments in hadronization induced asymmetries. 



beam energy. The combination of excellent vertex-region 
tracking resolution (typically 65 /xrad at 150 GeV/c) and 
excellent momentum resolution (cFp/p ~ 0.5% for a typ- 
ical 100 GeV/c track) gave a A c + mass resolution (a) 
of about 9 MeV/c 2 at all x F . Tracks above 22 GeV/c 
traversed a Ring-Imaging Cherenkov (RICH) and were 
identified as 7r, K, or p §. The RICH distinguished K± 
from 7T^~ up to 165 GeV/c. The proton identification effi- 
ciency was greater than 95 % above proton threshold (« 
90 GeV/c) and greater than 90 % above kaon threshold 
(w 45 GeV/c). 

An important innovation in SELEX was the use of an 
online topological trigger to identify charm. The hard- 
ware trigger was loose, requiring > 4 charged hadrons 
in the forward 150 mrad cone and > 2 hits from posi- 
tive track candidates in a hodoscope after M2. About 
1/3 of all inelastic interactions satisfied this trigger. The 
software trigger made a full vertex reconstruction of the 
beam track and all tracks in the M2 spectrometer (high 
momentum tracks) to test that hypothesis that they all 
came from a single primary vertex. A software-adjustable 
X 2 cut of 8.5 selected candidate events that were inconsis- 
tent with a single primary vertex, i.e., events that might 
have a downstream decay. This gave a rejection factor of 
about 8 in the data volume. Studies with the software 
trigger turned off show that this costs about a factor of 
2 in charm efficiency, for a net filter enhancement of 4. 

III. DATA ANALYSIS 



II. APPARATUS DESCRIPTION 

The SELEX spectrometer is a three-stage magnetic 
spectrometer designed to cover charm production in the 
forward hemisphere. The beam hadron was identified by 
a transition radiation detector as pion or not-pion and 
tracked in 8 beam silicon strip detectors before inter- 
acting in a set of 5 target foils, 2 copper and 3 carbon 
totalling 5% of a proton interaction length. Interaction 
products were tracked in a set of 20 vertex silicon detec- 
tors arranged in 4 sets of planes rotated by 45°. Single- 
track efficiencies for each system exceeded 98%. Data 
were taken with negative beam at 615 GeV/c (half 7r, 
half £) totalling 12.3 x 10 9 inelastic interactions and 
positive beam at 540 GeV/c (92% protons) totalling 2.8 
x 10 9 inelastic interactions. The layout of the spectrom- 
eter can be found elsewhere Q . 

Two Proportional Wire Chamber (PWC) magnetic 
spectrometer systems analyzed tracks with momenta 
above 3 GeV/c (Ml spectrometer, 12 planes) or 15 GeV/c 
(M2 spectrometer, 14 planes). The beam region of the 
PWCs was covered with double-sided silicon detectors to 
give high-resolution momentum measurement up to full 



A general data reduction pass found reconstructable 
charged tracks and identified the interacting beam track. 
Primary and secondary vertices were found by geomet- 
ric reconstruction without regard to particle identifica- 
tion. The high resolution of the beam and vertex sys- 
tems, along with the thin production targets, gave typ- 
ical primary vertex resolution of 270 /jm and secondary 
vertex resolution for these A+ events of 560 /j,m for a 
mean A+ momentum of 230 GeV/c. 

All the pKir candidates used in this analysis were se- 
lected by the following requirements: (1) primary and 
secondary vertex fits each had x 2 /dof < 5; (2) the 
primary-secondary vertex separation significance (L/cr) 
had to be greater than 8, where a is the quadrature sum 
of the primary and secondary vertex errors; (3) the to- 
tal momentum vector of the decay tracks had to point 
back to the primary vertex within errors; (4) (A c ) A+ 
events were required to have positive RICH identification 
for the (anti)proton and the (K + )K~. The pion could 
be any track; (5) secondary vertices which occured inside 
any target were removed; (6) each secondary track was 
extrapolated back to the primary vertex to determine its 
transverse miss distance. The second-largest transverse 
miss distance had to exceed 20 /xm. 
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Events were selected by the above cuts and sorted by 
the primary beam tag: E~ or n~ from the negative beam 
and protons from the positive beam. Figure |l| shows 
the full SELEX statistics for the A c — > pK + n~~ and 
A+ — > pK~n + mass plots where signals for all beams 
are combined. Gaussian fits are shown for illustration 



700 




2.2 



2.4 

Mass [GeV/c ! ] 



2.2 



2.4 

Mass [GeV/c 2 ] 



FIG. 1. Mass distributions for A c ~ (left) and A+ (right) 
where signals for all beams are combined. Gaussian fits are 
shown for illustration only. 

only. The background shapes are smooth on the low mass 
and high mass side of the signal regions. The signal yield 
for bayons and antibaryons are determined by a side- 
band subtraction method which assumes a linear mass 
resolution dependence for the momentum spectrum. The 
number of background events estimated outside the sig- 
nal window is subtracted from the number of events in 
the signal window. The signal window corresponds to a 
3.33-er mass width and the background is 10-a centered 
at 2.285 GeV/c 2 . We find 1979 ± 71 signal events for A+ 
and 520 ± 60 signal events for A c 



IV. ASYMMETRY ANALYSIS AND X F 
DISTRIBUTIONS 

The data yields were corrected in each xp bin by the 
total acceptance (geometrical acceptance and reconstruc- 
tion efficiency) . The simulation used Monte Carlo gener- 
ated A+ (or c.c.) events with phase space decay to the 
pKn final state. The simulated trajectories were con- 
verted to hits in all tracking detectors, including multiple 
Coulomb scattering but not secondary interaction losses. 
The hits were embedded in the hit banks of real data 
and all tracks were reconstructed. This overestimates 
slightly the tracking inefficiency, but this effect is small. 
The RICH simulation has been tested against a set of 
A — > p + ir~ events taken without a RICH cut. Agree- 
ment is excellent; the proton identification efficiency av- 
eraged over the momentum spectrum is 95%. 

The simulated events were analysed with the same 
tracking and reconstruction procedure used in the data. 
The whole set of cuts used to extract the signal was ap- 
plied here as well. The acceptance is the ratio between 
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0.15 - 


0.85 


2.65 ± 0.44 


0.15 - 0.675 


2.2 ±0.8 


£" 


0.125 - 


0.925 


2.45 ±0.18 


0.125 - 0.825 


6.8 ± 1.1 


P 


0.15 - 


0.85 


2.33 ±0.30 


0.15 - 0.60 


no fit 



TABLE I. Summary of measured n-value from the fit 
(1 — x F ) n for A* and A c . The errors are statistical only. 



the number of events reconstructed by the number of em- 
bedded events generated. SELEX acceptance does not 
depend on the transverse momentum of the charm par- 
ticle up to p\ = 7 (GeV/c) 2 . Therefore we consider only 
the acceptance as a function of xp. An important issue 
for this analysis is the relative efficiency for the baryon 
and antibaryon decays. As seen in figure |2|, there is no 
discernable difference in the acceptances. The Xp dis- 
tribution of the acceptance-corrected number of events 
from each beam is shown in figure |^. Because the beam 
flux is the same for baryon and antibaryon, these fig- 
ures compare the relative production cross sections for 
these states. The curves shown are fits of the standard 
parametrization (1 — xf) u to the data. The results are 
summarized in table |. The relative sizes of the differen- 
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FIG. 2. Acceptance of A c reconstruction. 

tial cross sections for baryon and anti-baryon production 
are clearly different for pion and baryon beams. The pion 
beam valence quark content is baryon/ antibaryon sym- 
metric. Relative cross sections for A+ and A c produc- 
tion from pions are comparable in shape and magnitude 
over the entire xp range. Such behaviour is consistent 
with the leading particle picture, since both charm states 
share a valence quark in common with the beam. For 
£~ and proton beams A+ production is much stronger 
than A c production, especially at high xp- In this case 
there are no valence antiquarks, and antibaryon produc- 
tion appears to be strongly disfavored. The four distri- 
butions that show leading behaviour are consistent with 
a common (1 — xp) n dependence having n ~ 2.5, much 
harder than the n for leading meson production from 
pion beams |^||. Note that the high-statistics £~ dis- 
tribution shows some interesting structures. Clearly A+ 
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production is favored over A c , and the enhancement in- 
creases with xf ■ The p 2 spectrum for A+ production by 
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FIG. 4. A c pt distribution of corrected number of events 
for 7T~ (left), proton (center) and S~ (right) beams. 



FIG. 3. A c xf distribution of corrected number of events 
for 7T~ (a), proton (b) and E~ (c) beams, (d) shows the A^ 
xf distribution for E~ production at low pf and high p\ . 

the three different beam hadrons is shown in Fig. |^. For 
pf < 2 (GeV/c) 2 all three beams are consistent with a 
smooth exponential spectrum e~ bPt with the same pa- 
rameter b « 1.1 (GeV/c) -2 . At larger p 2 values, seen 
only in the high-statistics S~ data, the exponential gives 
way to a power-law behavior as expected in QCD and as 
seen m D ± '° production g. There is no indication in the 
large- x f events of an enhancement at small p 2 as sug- 
gested by the intrinsic charm picture. However, when we 
look at the xf behaviour of events after the slope change 
at p 2 > 2 (GeV/c) 2 , we see a significant difference in the 
xf behaviour. The events with large p 2 have a softer xf 
spectrum than those at small p 2 , as seen in figure || (d). 
The hadroproduction asymmetry A is defined as: 



A = 



(1) 



where a c (ctc) is the production cross-section for the 
charm particle (anti-charm particle) in study. 

We calculated the values of the asymmetry and the 
errors using the maximum likelihood method with no 
limit at A = ±1. Events in each xf bin follow the Pois- 
son distribution. Therefore, the likelihood function is a 
product between two Poisson distributions, one for the 
background region and the other for the signal region. 

The results are presented for the three beam hadrons 
in figure]^. As is already clear from the xf distributions, 
A c production in £~ and proton beams shows a strong 
asymmetry favoring baryons over antibaryons at all Xp ■ 



The proton asymmetry rises sharply toward unity already 
at small x p , reflective of the paucity of A c production by 
protons. The S~ asymmetry has a smoother approach to 
the limit at +1. In the case of the ir~ beam we see only a 
slightly positive and constant asymmetry throughout the 
xf range. Our result agrees well in the overlap region 
with recent E791 data l[10|. 



V. SYSTEMATIC ERRORS 

Systematic errors in kinematic analysis may arise from 
several sources. We list here the ones that have been 
considered. In all cases systematic effects arc found to 
be small compared to statistical errors. 

• xf binning effects 
The xp distribution is relatively smooth, but there 
are structures in the high statistics data. We have 
shifted the bin limits for all xp distributions by 
0.025 (half a bin) and refit the data. No effect 
in the character of the distributions is observed. 
For example, all the structures seen in Fig. 3(c) are 
also present in the shifted distribution. Changes 
in the n- value are consistent with purely-statistical 
effects. Therefore, we quote a conservative An for 
the fit of 0.1, even though no systematic effect has 
been seen. 



• Mass window size 
The signal window in the mass distribution is taken 
as ±30 MeV/c 2 around 2.285 GeV/c 2 . We have 
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FIG. 5. Asymmetry for A c production by 7r (top), S 
(center) and proton (bottom) beam. 

changed the window up to ±35 McV/c 2 and down 
to ±25 MeV/c 2 . The systematic error for changing 
the mass window compared to statistical error is 
negligible. 

• Acceptance corretion 

The acceptance shown in Fig. 2 is computed using 
the method of adding hits from simulated charm 
decay tracks into the data banks from an observed 
event. The modified event is rcanalyscd to sec if 
the charm state is reconstructed. There could be a 
shift in the corrected xp distribution if the accep- 
tance calculation does not represent truly the xp 
variation of the apparatus efficiency. In E781 the 
dominant efficiency variation at low xf is purely 
geometric: both proton and kaon pass through the 
fiducial volume of the RICH counter. A second 
geometric effect is the related one: the charm par- 
ticle live long enough to pass the vertex significance 
cut. The tracking characteristics of the spectrome- 
ter match those predicted from the simulation very 
well over all momenta studied here. The K s mass 
is the same for all observed momenta, both in data 
and simulation. We believe that there is no system- 
atic bias from the simulation and assign no error. 

• Trigger Bias in Asymmetry 

The XF-dependent effects tend to cancel in the 
asymmetry measurement. However, a new system- 
atic effect must be considered: Trigger bias. The 
E781 trigger required counts on the postive-bend 
side of a trigger hodoscope after the second bend- 
ing magnet. This might favor A+ production over 



At- Wc looked at the asymmetry of K s produc- 
tion, for times when the tt + is found in the high- 
momentum spectrometer compare to the converse. 
We see a maximum 3% asymmetry. We also look 
at the charm baryon data itself. If we exclude the 
charm tracks, then in 97% of the triggers there are 
still enough hits in the hodoscope to trigger the 
event without the charm decay tracks. We conser- 
vatively assign a 3% systematic error to the asym- 
metry from trigger bias. 

VI. SUMMARY 

The SELEX experiment reconstructed about 2500 
A c — > pKir from 600 GeV/c X - and ir~ beams and 
a 540 GeV/c proton beam. We observed that the ip- 
dependence of A+ production is similar for all three 
beams, much stiffer than for meson production. Both 
baryon beams show a strong enhancement of the pro- 
duction of charm baryons over antibaryons, while the two 
are produced comparably from a pion beam. These re- 
sults are in general agreement with color-drag models of 
hadronization. The differences in behaviour for protons 
and S~ illustrate significant non-perturbative effects in 
production dynamics at these energies. 
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